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Abstract 

^ , We have investigated the energy loss of hot electrons in metallic graphene by means of GHz 

noise thermometry at liquid helium temperature. We observe the electronic temperature T (x V 
at low bias in agreement with the heat diffusion to the leads described by the Wiedemann-Franz 
law. We report on T oc yV behavior at high bias, which corresponds to a T'^ dependence of the 
cooling power. This is the signature of a 2D acoustic phonon cooling mechanism. From a heat 
equation analysis of the two regimes we extract accurate values of the electron-acoustic phonon 



^]^ ■ coupling constant S in monolayer graphene. Our measurements point to an important effect of 

cn 

^vq I a strong and firm support to the rising field of graphene bolometric detectors. 
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lattice disorder in the reduction of S, not yet considered by theory. Moreover, our study provides 
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Many remarkable properties of graphene stem from the chiral Dirac-fermion nature of 



QE 



carriers and the associated semi-metaUic density of states [1|, |2| . Graphene is also peculiar 
from the phonon point of view with a large optical phonon (OP) energy of 160-200 meV 
and a weak coupling of electrons to acoustic phonons (APs) 3|-[6j. The later is due to the 
non-ionic character of the lattice and the 2D nature of the phonons. Both effects conspire 
to the weak electron-phonon scattering responsible for the high intrinsic carrier mobility of 
graphene up to room temperature [7|. The OP-scattering reappears however in the high- 
ield transport, where OP emission by energetic electrons gives rise to current saturation 
8|, |9|. By contrast, extrinsic disorder effects usually obscure AP-scattering which remains 
thus more elusive and still debated [10|. The AP-coupling becomes however prominent in 
heat transport as it constitutes the bottleneck in the cooling process of hot electrons to the 
substrate. Heat transport experiments in graphene offer thus the possibility to characterize 
electron-AP coupling in detail. On the application side, the weak AP coupling opens an 
avenue for the realization of fast and sensitive bolometric detectors working in a broad 



frequency spectrum ranging from microwave to optics lll-ll5l|. Also in this respect an in- 
depth and direct characterization of hot-electron effects will prove useful. 

In this letter we report on the interaction between electrons and acoustic phonons in 
mono-layer graphene. To this end, we have investigated the bias dependence of the electron 
temperature deduced from radio frequency (RF) shot noise measurements. Our experiment 
provides a direct access to the electron-phonon coupling constant. We rely on the relation 
Si = AksTe/R between the current noise spectrum S"/, the average electronic temperature 
Tg, and the sample resistance R. Johnson-noise thermometry has proven useful to study 
carbon nanotubes ly, [l7| and more recently bilayer-graphene [l8|. From an experimental 
point of view, we measure noise in the MHz to GHz band to overcome the environmental 
1// contribution while keeping quasi-stationary response conditions [l9(]. At sufficiently 
high bias, we find Tg oc ^A^, where V is the applied voltage. This is a signature of 2D 
acoustic phonons. We observe it in diffusive samples where the I — V characteristics are 
linear and optical phonon scattering is absent. At low bias we recover the T^ oc V behavior, 
expected for heat conduction to the contacts, and described by the Wiedemann-Franz law. 
Both regimes are well explained by the heat equation. By fitting it to our measurements 
we extract the coupling constant S as function of carrier density for samples on different 
substrates. 



sample 


L xW {jj.m^) 


R(k^ ) 


BNl 


2.2 X 5.7 


2.8-3.8 


BN2 


2.2 X 2.7 


1.3-2.3 


CVDl 


1 X 1 


1.67 



TABLE I: Characteristics of the graphene samples. L is the sample length, W the sample width, 
and R the drain-source resistance. 

We present data of three representative devices: two large exfoliated graphene flakes 
deposited on hexagonal boron nitride (hBN) and a smaller sample on Si02 made from 
chemical vapor deposition (CVD) graphene. The samples, labelled BNl, BN2, and CVDl 
respectively are sketched in Fig{Tl Sample BNl and BN2 consist of a stack of a monolayer 
graphene and a 30 nm-thick hBN platelet on a Si/Si02 substrate. The hBN platelet of 
approximately 10 /im diameter was produced by exfoliation of a high quality hBN powder 
(St. Gobain "Tres BN"). The graphene flake was also obtained by exfoliation and placed on 
top of the hBN using a wet transfer techniqu e |20i| . CVDl was made from a graphene sheet 
grown on copper, transferred to the Si/Si02 [21|. The samples were produced by means of 
e-beam lithography and dry etching. Doped silicon was used as a back gate, with a 1 fun 
thick oxide, and each sample was embedded in a coplanar wave-guide adapted for GHz 
frequencies 22|. Palladium metallizations ensure low ohmic contacts to graphene. Tabled 
lists the samples' dimensions and electrical characteristics. 

Measurements were carried out in liquid helium at 4.2 K to ensure a cold phonon bath. 
The ampliflcation line, of bandwidth 0.8 GHz, is similar to the one in Ref. [16]. It in- 
cludes a cryogenic low noise amplifler, calibrated against the poissonian noise of a 460 Ohms 
Al/AlOx/Al tunnel junction. The samples' I-V characteristics were measured via the volt- 
age drop across the bias resistance. In situ current annealing was performed to improve 
the samples' electrical characteristics. The carrier density n^ of samples BNl and BN2 are 
deduced from the gate dependence of the DC conductivity in a range of V^ = ±18 V(BNl) 
and Vg = ±28 V(BN2). Sample BNl stayed p-doped over the whole gate range whereas 
BN2 could be swept through the neutrality point. BN2's mobility is on the order of 
3000 cm^V^^s^^, about a factor 10 larger than that of BNl. The conductance of CVDl 
was independent of Vg which indicates a highly doped regime; we estimate the carrier con- 
centration on the order of a few 10^^ cni^^ from Hall-bar measurements on similar CVD 



sheets. 

We turn now to the presentation of measurements. Fig|2]^a) shows typical current noise 
spectra for different bias voltages. At low bias the spectrum is flat, whereas at higher bias a 
1// contribution clearly arises. Our broad-band setup allows us to quantitatively separate 
the white shot noise contribution Sj from 1// noise by fitting the spectra with a S"/ + C/ f 
law. We find C oc /^ in accordance with the Hoodge law 23|. This procedure is especially 
important for small samples where the 1// contribution is larger. Figs|2]^b) and (c) show 
the current I and the shot noise Sj of sample BNl as a function of bias for three gate 
voltages. We restrict our study to biases where I-V characteristics are linear as the lack 
of current saturation is an indication of the absence of optical phonon scattering 81]. This 
is consistent with a transport scattering length (< 10 nm) much smaller than the optical 
phonon scattering length (> 200 nm) [S\. The shot noise is strongly sub-linear with respect 
to bias voltage and depends on gate voltage, with larger noise at higher carrier concentration. 
We extract the bias dependence of the average electronic temperature Te(y) = SiR/Aks 
from the measured shot noise. The results for sample BNl at five gate voltages are displayed 
in Figj3l The electronic temperature increases with bias voltage and reaches several hundred 
Kelvin at high bias. This value is well above bath temperature but still below the optical 
phonon energy of ~ 200 meV. We find that Te is weakly dependent on gate voltage and 
the different curves closely follow the Tg oc \/V law (dashed line) predicted for 2D acoustic 
phonons. Deviations are however observed at low bias where data rather follows a T^ cc V 
law. At the low end of the investigated bias range, the quasi-equilibrium approximation 
required to define an electron temperature may fail due to large electron-electron scattering 
lengths 2J|. From a Fano factor analysis, we estimate this mesoscopic regime to arise below 
a few tens of mV |33| . 

The electron temperature is determined by the balance between the Joule heating and 
the cooling powers at play 251]. Joule heating per unit area is P = V'^/RWL, while cooling 
is provided by two mechanisms: The first is electron heat diffusion to the leads, the second 
electron- AP interactions. At sufficiently high bias the diffusive contribution can be neglected, 
i.e. the cooling power is given by T,{T^ — T^f^), where Tp^ is the phonon temperature. In 
graphene 5 = 4 has been proposed to account for the purely two dimensional APs [4, l6| . In 
general 6 varies from 3 to 6 26|]; the lower value 5 = 3 is observed in ID systems such as 



carbon nanotubes p/7[ and 5 = 5 is typically used for 3D metals [24 



271 ]. For a quantitative 



analysis of our data, we solve the heat equation: 



"^"^'^'^'^"^ -^ + WS[T^(x)-Ti] , (1) 



2R dx^ R 

where C = 7r^A;^/3e^ is the Lorenz number and x denotes the coordinate along the graphene 
channel. We obtain an analytical solution in the case of cold contacts and cold phonon 
bath [Tph = 0). The solution depends on a single free parameter S, which sets the two 
characteristic scales of the system: the temperature Ts = {V'^/RLWT.y^'^ and the voltage 



Vs = C/y/ARLWTi. Ts = (P/S)^/^ is the maximum temperature reached in the absence 
of electron heat conduction, whereas Ve defines the crossover voltage between the electron 
cooling at low voltages and phonon cooling regime at high bias. The temperature profile is 



pseudo parabolic at low bias [25| and evolves toward a uniform temperature T{x) = T^ at 
high bias. The spatial average Te =< T{x) > is then calculated numerically. 

In FigJH we plot the temperature data in the form T^/P as a function of V with the 
corresponding error bars. Due to a more pronounced 1// contribution, the error increases 
with increasing bias. The T^/P representation allows for a better comparison with theory 
and puts the emphasis on the density dependence of S. The plateaus at high bias and the 
dips at low bias reflect the Tg oc \/V and Te oc V regimes observed in the shot noise data 
(FigJ3]). The measurements are well fitted by the solution of Eq.([T]) taking S as the only 
free parameter (solid lines). We find S < 2 mW/m^/K^ and a clear dependence on carrier 
density in BN samples, and S ^ 0.42 mW/m^/K^ in CVDl. Since the dip at low bias marks 
the effect of electron heat diffusion to the leads, the crossover is more pronounced in the 
shorter CVDl sample. This can be seen in FigJl](a) where the width of the dip is larger 
and agreement between data and theory is more quantitative. The fact that theory, with a 
single free parameter, accounts for both the plateau and the dip is a strong confirmation of 
the model and the underlying 2D acoustic phonon mechanism. FigJl](c) shows the carrier- 
density dependence S(ns) and emphasizes its sensitivity to carrier mobility and disorder. 
This observation is consistent with the smaller values of S in the CVDl sample. 

The afore mentioned STg^ law for electron-AP cooling is similar to a Ei<-TA law for 
phonon black body radiation to the substrate or the bath (Tq). This mechanism would 
occur in a hot phonon regime where Tph — Te ^ Tq. However the corresponding coupling 
constant T,k is about 3-4 orders of magnitude larger than S measured in our experiments 



28 



29| . Both mechanisms appear in series in the electron cooling process, thus giving rise 



to an effective coupling E = SS/^/(S + Hk) — S. This justifies our earlier hypothesis 
TA ^ T^ in Eq.([T]). In order to confirm our cold acoustic phonon approximation we have 
performed Raman spectroscopy measurements in a similar sample. As already reported 



by several groups 30|, |3l|, we find that the position of the G (resp. 2D) band shows a 
downshift of 0.024 cm~^/K (resp. 0.051 cm^^/K) when increasing the cryostat temperature 
(FigJSt^a)). When the sample is biased, a similar downshift is observed, which is attributed 
to an increase of the AP bath temperature (Figl5](b)). This behavior was confirmed both 
at room temperature and at 100 K. We estimate that the corresponding increase of the AP 
bath temperature remains below 30 K/V in all cases. As compared to electronic temperature 
elevations (Figj3]), these values entail negligible corrections into Eq.(Il]). 
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We now come to the discussion of the results. According to theory |4, |6[, longitudinal 
acoustic phonons (LA) are coupled to electrons via the deformation potential D (10-30 eV). 
Theory predicts a cooling power P = S/^^T^ with T^la oc D'^^fnl in the non degenerate 
metallic regime and a different PiT^ law at high temperature. Taking D = 10 eV, this gives 



Y^LA — lOyn^/To^ mW/m^/K^. Our experiment confirms the P = ST^ dependence in the 
metallic regime, but we find a smaller coupling constant, S = 0.5 mW/m^/K^ for BNl and 
E = 2 mW/m^/K^ for BN2 at Ug = 10^^ cm~^. The discrepancy in S cannot be explained 
by experimental uncertainties. Our experiment suggest the effect of lattice disorder, which 
is the limiting factor for the mobility, as a possible cause of the reduction of the electron- 
phonon coupling. As a matter of fact, lattice disorder, which is known to affect phonon 

n ri 

lifetime, is not taken into account by theory for the metallic regime |^ |6|]. In the high- 
temperature regime however, a recent theory shows that disorder mediated phonon-phonon 



collisions may be important 32 1. 



In conclusion, we have studied the interaction between electrons and acoustic phonons in 
diffusive graphene by measuring the energy-loss of hot-electrons. We find that the cooling 
power due to acoustic phonons follows a ST^ law. We find that S increases with carrier 
concentration, but the value, S ^ 2 mW/m^/K'^, is approximately one order of magnitude 
smaller than predictions for LA-phonons. Our work is a motivation for further theoretical 
investigation that takes into account the role of lattice disorder on phonons. Beside its 
implication for electron-phonon physics, our study of thermal noise is of direct relevance for 
the performance of graphene detectors. 
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FIG. 1: Shot noise measurement scheme (color onhne) with a sketch of the sample the electrical 
environment. The samples and the cryogenic low-noise amplifier are cooled at 4.2K. The SEM 
pictures refer to samples BNl and CVDl: BNl consists of an exfoliated graphene flake (red), 
deposited on an exfoliated h-BN platelet (green), itself deposited on a Si/Si02 used as a back- 
gate. CVDl was fabricated from CVD grown graphene deposited on Si/Si02. Both samples are 
contacted with Pd electrodes (yellow). 
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FIG. 2: (a) Typical excess noise SY^{V) spectra in sample BNl; it is a white noise with a superim- 
posed 1// contribution and fitted by SJ^ = S/ + C/f laws (solid lines), (b) and (c) show the I{V) 
and Sj{V) data for different gate voltages in sample BNl from which we deduce Tg. The circles 
point out the I{V) and Si{V) values of the spectra shown in (a). 
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FIG. 3: Electronic temperature in sample BNl as function of voltage bias for a set of gate voltages. 
Te = SjV/{4:kBl) is deduced from the Si{V) and I{V) data shown in Figj2]^c) and (b). Unlike 
Sj{V), Te{V) is nearly independent of gate voltage and closely follows the Tg oc VV law (black 
dashed line) expected for 2D phonons. Deviations are observed at low bias where aT^ ccV behavior 
is found. 
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FIG. 4: (a) and (b): Electron temperature of sample CVDl (a) and BNl (b) plotted as T^{V)/P, 
where P is the Joule heating per unit area P = V'^/RLW. The plateau at high bias is at a value 
T^/P ~ 1/S. The dip at low V is due to electron heat diffusion to the leads. Dashed lines are 
one-parameter fits with S as free parameter, (c): S as function of carrier density Ug for samples 
BNl and BN2. 
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FIG. 5: (a): Raman shift of the 2D-band of a graphene sample similar to CVDl, as function of 
the cryostat temperature Tq. The slope is (—0.051 it 0.008) cm^^K^^. Excitation laser at 532 nm 
, P = 25 kWcm^^. No laser power dependence was observed in this range, (b): Raman shift of 
the 2D band as function of the bias voltage for Tq = 100 K (blue symbols) and Tq = 300 K (red 
symbols). The slopes are (—0.3 it 0.3) cm~^V^^ and (—0.9 it 0.3) cm^^V^^ respectively. The bias- 
induced phonon heating is therefore below 30 K/V. Similar results and conclusions were drawn 
from measurements on the G band (not shown). 
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